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Abstract. The notion of the Radon transform Rh on the Heisen- 
berg group ]H[„ was introduced by R. Strichartz and inspired by 
D. GeUer and E.M. Stein's related work. The transversal Radon 
transform Rt integrates functions on R™, m > 2, over hyperplanes 
meeting the last coordinate axis. In fact, Rh is a particular case 
of Rt, corresponding to m = 2n + 1. We obtain new boundedness 
results and explicit inversion formulas for both transforms on 
functions in the full range of the parameter p. We also show that 
Rh and Rt are isomorphisms of the corresponding Semyanistyi- 
Lizorkin spaces of smooth functions. In the framework of these 
spaces we obtain inversion formulas, which are pointwise analogues 
of the corresponding formulas by R. Strichartz. 



1. Introduction 

Let M.n = C" X M be the Heisenberg group with the multiphcation 
law 

{z, t) o (C, r) = (2; + C, t + r - - Im{z ■ ()) 

and let 

(LI) {RHf){z,t)= [ /((z,t)o(C,0))dC= / f{C,t-llm{z-0)dC 

be the Radon-like transform, which was introduced by Strichartz in his 
remarkable paper [Str2l| . He called Rh the Heisenberg Radon transform 
and noted that definition (11.11) stems from the previous work by Geller 
and Stein [GSlt IGS2] . The Heisenberg Radon transform represents an 
interesting object from the point of view of both harmonic analysis and 
integral geometry. Related transforms in the more general context of 
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nilpotent Lie groups and Siegel domains were studied by Felix |Felt 
He and Liu [HLTl [HL2] . Peng and Zhang [PZ]; see, also [HiH 
EiaiLHllNT]. 

A variety of deep results in the area of harmonic analysis on the 
Heisenberg group can be found in fundamental works by M.L. Agra- 
novsky, C. A. Berenstein, D.C Chang, M. Cowling, G.B. Folland, D. 
Geller, A. Koranyi, D. Miiller, D.H. Phong, F. Ricci, E.M. Stein, S. 
Thangavelu, J. Tie, and many other authors; see, e.g., [BCTt [UK| IKRt 
IKorl IMRSt ISt2l IThan] , and references therein. 

Owing to |Str2l p. 386, 387], in a slightly different notation we have 
the mixed norm estimate 



l<P<Po, Po = l+l/i'2n+l), r = p/{2n+l-2np), l/p+l/p' = l. 
The paper [Str2j (see pp. 386, 399) also contains the inversion formulas 



(L3) R-^' = (4vr)-^"(^)"i?H (^)", Rh' = ^-'J^2, 



where JF2 denotes the Fourier transform in the t-variable, is the 
Fourier transform in all variables, and {Jf){x,y) = f{—2y~^x,y). 

Our research is motivated by the following questions, that stem from 
[Str2] and remained open in afore-mentioned publications. 

Question 1. Is RhJ well defined a.e. for / G L^(EI„), when p > p^l 
If yes, then what can we say about the boundedness oi Rh 'i^ this case? 

Question 2. What is the substitute of ([L3D for / G ^^(e^)? 

Question 3. To what class of functions are formulas flL3l) applicable 
pointwise? 

In the present article we answer these questions in a more general 
context of the so-called transversal Radon transform Rt-, which in- 
cludes Rh as a particular case. Our method differs in principle from 
those in the previous related works j5tr2l [Fell [Fe2l [HLTl [HL21IPZ] . To 
be more specific, let us rewrite (11. ip by setting z = u + iv X = i + 'i'Vi 
where u, v,^,!] G M". We obtain 



Then we identify /(^ + ir],T) with a function f{x) on M™, m = 2n+l, 
where 



(1.4) x = {x',Xm)eW, x' = {^,r])eW-\ x^ = rGR, 



(1.2) 
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and set 

(1.5) a = ^{-v, u) e b = teR. 
This gives 

(1.6) {RHf){u + iv,t) = {RTf){a,b), 
where 

(1.7) (R^f)(a,b)= [ f{x',a- x' + b)dx'. 



Integral (11.71) can be regarded as a Radon transform of a function / on 
associated with the hyperplane 

(1.8) h = {x = {x, x„) G R"" : Xm = a-x' + b}, 

which is transversal to the last variable. Following Strichartz [Str2, 
p. 385], we call (11.71) the transversal Radon transform. It resembles a 
special case of the affine Radon transform studied in Gelfand's school 
|GGGl IGGV] and the parametric Radon transform in ^Ehrj . To distin- 
guish notation, we write 

(a, 6) G M*" = M""^ X M. 

Note that m in (ll.7p may be an arbitrary integer > 2. The case of 
m odd agrees with the Heisenberg Radon transform Rh- Thus, basic 
results for Rjj will follow from those for Rrp. 
Main results. We invoke the classical Radon transform 

(1.9) {Rfm= [ f{x)dh{x), 

Jh 

where integration over the hyperplane h in M"* is performed against 
the usual Lebesgue measure. Suppose that h is parameterized as 

(1.10) /i = {x G : x-e = t}, {6, t) G = S""-^ X M, 
^m-i ]3gi]2g tihe unit sphere in M"^. Then 



(1.11) iRf)ih)^iRf)ie,t)= fiy + te)dy. 
The corresponding dual transform is 

(1.12) {R*tp){x)= I ip{e,x-e)de. 

The duality relation 

{Rf){e,t)i){e,t)dedt= [ f{x){R*i)){x)dx 
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holds provided that at least one of these integrals is finite when / and 
ijj are replaced by |/| and IV'I, respectively. 

Integral transforms (II. 9p and fll.l2p . are the most familiar to the 
reader. They were introduced by Radon |Rad] and studied in numer- 
ous publications together with their fc-plane generalizations; see, e.g., 
[BiH iDiil [GGGl EU Kitl [Mill Eil [Ml [QCKl [RKl [RuTl [SSW] and 
references therein. 

Unlike fll.lOp . that parameterizes all hyperplanes in M™, f ll.Sp ex- 
cludes hyperplanes parallel to the If both parameterizations 
are available, then the standard Calculus yields 

(1.13) {Rf){h) = ^lTW{RTf){aM- 

Equation fll.l3p means that basic properties of Rh and Rt can be 
derived from known facts for R. However, as indicated by Strichartz 
|Str2l p. 385], there are essential distinctions between Rt and R. This 
remark inspires us to give direct proofs, that might be instructive. 

Theorem 1.1. Let Rea > 0, 

~r((a+m-l)/2)- 
The following equalities hold: 

(1-14) / (^T/)(a,&) I ,^^„,^,, =A / f{x)\xrUx, 
(1.15) 

fj, m^'-'dadb f , \x\^-' , 

provided that either side of the corresponding equality is finite when f 
is replaced by \f\. 

Equalities (11.140 and (11.151) are proved in Section 2 and play a key 
role. As we shall see in Sections 3 and 4, (I1.14p paves the way to 
the "right" definition of the dual transversal Radon transform, the cor- 
responding analytic families of fractional integrals of the Semyanistyi 
type, and a variety of explicit inversion formulas for Rt. Equality 
(I1.15P answers the question about existence of the transversal Radon 
transform of functions and provides some boundedness results. In 
particular, by Holder's inequality, (11.150 implies the following corollary 
(cf. [Rull Corollary 2.4]) . 

Corollary 1.2. If f e Lp(M"^), 1 < p < m/(m - 1), then {RTf){a,b) 
is finite for almost all (a, b) G M™". Moreover, for any a > 1 — m/p' , 
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l/p' = 1, we have 

\bY^~^dadh 

liRrfKa, b)\ ^^^ ^^^^ ^ 52)(a+m-i)/2 ^ ^\\J\m . 
where c = c{a,p,m) = const < oo. 

Corollary 11.21 leads to the following result for the Heisenberg Radon 
transform (just set m = 2n + 1). 

Theorem 1.3. Let = 1 + l/2n. If f e Lp(H„), I < p < pi, 
then {RHf){z,t) is finite for almost all {z,t) G ]HI„. Moreover, for any 
a > 1 - {2n + l)/p', 

Itl^-^dzdt ^ 

|Lf{H„)- 





t 


"^-^dzdt 


(1 + 


z 


2 _|_ ^2^(a+2n)/2 



Clearly, px = l+l/2n in Theorem II .31 is greater than pq = 1+1/ (2?2+l) 
in fll.2l) . This gives the answer to Question 1. Note also that the 
condition p < m/{m — 1) in Corollary 11.21 (resp. p < pi in Theorem 
11.31) is sharp. This follows from (11.131) and the well-known fact that if 
p > m/ (m — 1) and 

fix) = (2 + |x|)-'"/^(log(2 + \x\))-' (G LP), 

then {Rf){h) = oo; see, e.g., jSol IRu7] . 

Theorem 1.4. For m > 2 an a priori inequality 



\{RTf){aMdh 



holds if and only if 

(1.16) 1 < p < m/{m — 1), q = p', l/r = l — m/p'. 

This statement will be proved in Section 8. For m = 2ra + 1, it 
extends Strichartz's result (II. 2p to all 1 < p < pi, where the bound 
Pi = l + l/2n>po is best possible. 

A variety of inversion formulas for the transversal Radon transform 
and the Heisenberg Radon transform are presented by Theorems 14. 2^ 
SSI rm and Ol We develop severals methods involv- 

ing hypersingular integrals, powers of minus-Laplacian, and wavelet 
transforms. We also present the convolution backprojection inversion 
method for Rrf, f G Lp(M™), introduce the relevant ridgelet-like trans- 
forms and give examples of admissible wavelet functions. All inversion 
formulas for L^-functions are obtained in the maximal range of the 
parameter p. 



6 



B. RUBIN 



Another series of pointwise inversion formulas for Rt and Rh is 
obtained in the framework of the spaces of the Semyanistyi-Lizorkin 
type. For example, we introduce the space $(EI„), which consists of 
Schwartz functions (p on EI„ satisfying 

[ (p{C,T)T^dT = WkeZ+, VCgC". 
Jr 

We show that the Heisenberg Radon transform Rh is an automor- 
phism of $(]HI„) and obtain pointwise analogues of formulas fll.3l) . In 
particular, by Theorem 17.51 a function / G $(IHI„) can be reconstructed 
from (f = Rnf by the formula 

(1.17) /(C,r) = (-ir (47r)-2"(5"^M»(C,r), 

which differs from the first equality in (11.31) by the factor (—1)". To 
obtain these results we first establish more general relations for the 
transversal Radon transforms, including Riesz potentials. 

Some historical comments are in order. In 1960, Semyanistyi ^SeT| 
came up with an interesting idea to treat the hyperplane Radon trans- 
form and its dual as members of suitable analytic families of operators 

i?" and a G C, so that i?"/U=o = Rf and R"ip\o,=o = R*^- 
Semyanistyi's operators combine properties of Radon transforms and 
fractional integrals. They yield a variety of inversion formulas 

cj = (-A)("+"-i)/2 k'^iRf), 

where c„ = const and (— A)*^"+"^^^/^ is a power of minus-Laplacian, 
that can be realized in several ways. This philosophy was extended in 
different directions and proved to be useful in applications to convex 
geometry; see [nRl[Ru2] . [Ril4] -[Ril8]. \RZ \. 

In the present article we extend Semyanistyi's idea to the transversal 
Radon transform and demonstrate its benefits. 

We expect that some ideas of our article are applicable to Radon 
transforms on more general two-step nilpotent groups, as in |Felt IFe2[ 
\PZ\ IHLH IHL2] . We plan to address this topic in forthcoming publica- 
tions. 

Plan of the paper. Section 2 contains preliminaries and the proof of 
Theorem 11.11 In Section 3 we introduce analytic families of fractional 
integrals of the Semyanistyi type, associated to the transversal Radon 
transform Rt, and define the dual transversal Radon transform. This 
section forms a foundation for different inversion methods developed in 
Section 4. In Section 5 we introduce and investigate a function space 
$(M™) of the Semyanistyi-Lizorkin type associated to the transversal 
Radon transform. Section 6 is devoted to the detailed study of the 
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connection between the classical hyperplane Radon transform (11. 9p . 
the transversal Radon transform (II. 7p . the Heisenberg Radon transform 
(11.11) . and their duals. Results of this section are applied in Sections 7 to 
investigation of isomorphism properties and inversion of the Heisenberg 
Radon transform. Section 8 deals with mixed norm estimates of the 
transversal and Heisenberg Radon transforms. 

Acknowledgements. The work was started when I was visiting the 
Capital Normal University (Beijing) and the Guangzhou University in 
Summer 2009. I am grateful to Professors Zhongkai Li and Jianxun He 
for the hospitality, support, and discussions. They participated in some 
calculations related to Sections 5-8, and their results will be published 
elsewhere. I am especially thankful to Professor Dr. Rainer Felix for 
sending me his interesting paper |Fe2j . 

2. Preliminaries. Proof of Theorem 11.11 

2.1. Notation and auxiliary statements. We denote by ei, . . . , 6,^ 

the coordinate unit vectors in M™; 6 = (9,,...,9m) G S""-^; 9' 



m 

^ Qm-1. f)/ 

I, . . .9„ 



SI'-' = {9e : ft„ > 0}; 

am-i = 27r™/^/r(m/2) is the area of S*™""^; d9 stands for the usual 
Lebesgue measure on S"^~'. We recall the well-known Catalan formula 
from Calculus 



(2.1) 



/ f{9 ■a)d9 = am-2 f f{t){l - t2)(™-3)/2 



which holds for every a G S^^' provided that at least one of the 
integrals in (12.11) is finite when / is replaced by |/|; see, e.g., |ABR[ 
p. 216]. The notation C, C°°, and U for function spaces on R"* is 
classical; Cq = {/ e C{W^) : lim f{x) = 0}; 5(M'") is the Schwartz 

space of rapidly decreasing C°°-functions with standard topology. All 
function spaces on the Heisenberg group ]HI„ = C" x R are identified 
with the corresponding spaces on R^"+^. 

The Riesz potential J"/ on R™' is defined by 

(9 9^ (T-f\(.\ ^ / ^(^) ^ 2-7r™/2r(a/2) 

(2.2) (7 /)(x) = / — — -r— , 7m(a) - 



7m(a) Jm™ k - 1/1""-"' r((m-a)/2)' 

i?e a > 0, a — m 7^ 0, 2, 4, ... ; 

see |Stl] . The letter c stands for a constant that can be different at each 
occurrence; = {0,1,2,...}; denotes the relevant set of multi- 
indices. Given a real-valued expression A, let (A)+ = A'^ if ^4 > and 
{A)l = if A < 0. 
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Lemma 2.1. 

(i) Iffe L^S""-^), then 

(2.3) [ f{9)de= [ f{a)da, 



(ii) Conversely, if f & L^{W^~^), then 
(2-4) / na)cia=[ s{^)Z 



1 V6'm/ 



Proof. We write 

/ f{6)d9= / sin™^^ C(j (ic(j / f {a sin cu + Cm cos cu) da 

J s^-^ Jo J s^-'^ 

and then set s = tantu. This gives 

T -I- \ / f^n — \ 1 

da. 



(i + .2)W2y^, 



which coincides with (12.31) . Similarly, 

[ f{a)da= I s'^-^ds I f{sa)da 
Jr'"-i Jo J s-^-^ 

Jo COS^U; \ COSiU J Jg^-^ \6m^ dm 

□ 

We will be dealing with Riemann-Liouville fractional integrals 



(2.5) (/o+^)(r) = f^J^ {r-sr-'^{s)ds. 

The following fact from |Ru3] will be needed, where [a] denotes the 
integer part of a. 

Lemma 2.2. Let a > 0, 

POO 

/ s^ijj{s) ds = for all j = 0, 1, . . . , [a], 
Jo 

s'^\ip{s) \ ds < oo for some (3 > a. 



Then 



i(/o:»(^)i 



0(s") ifO<s<l, 
O(s^), 7 = a -min(l + [«],/?)< 0, if s > 1. 
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Furthermore, 



Jo 



s 



a+l roo 



poo 

/ 3°" ipi^s) log sds ifaeN. 
Jo 



More information about Riesz potentials and fractional integrals can 
be found in [Rul], [SKM]. 

2.2. Proof of Theorem II. IL We denote by / the left-hand side of 
f im]) . Then 

"^'^ ^ Ibl'^-Ub [ f{x',a-x' + b)dx' 



-1 (1 + |a|2)("+"^-i)/2 



m— 1 





a ■ 


x' 




o-l 


(1 + 


a 


2ya+m-l)/2 



By (12. 4p and (12.11) . the inner integral can be represented as 





a ■ 


X 




(1 + 


a 


2 


+ 


a ■ x' — 




2yQ+m-l)/2 



' ' 2 Js^-i r((a + m- l)/2) ' ' 

and (fTTil) follows. 

The proof of (11.151) uses the same idea. Denoting by / the left-hand 
side, we get 

f{x) dx 

By (12. 4p and (12.11) . the inner integral becomes 

io (l + r2t2)(a+— 1)/2^^' 

where r = The last integral can be computed by changing variables 
and using |GRy[ formula 3.238(3)]. □ 

3. Analytic families of fractional integrals and the dual 

TRANSVERSAL RADON TRANSFORM 

Below we introduce analytic families of fractional integrals associated 
to the transversal Radon transform. We start with (I1.14p and replace 
/ by the shifted function fx{y) = f{x + y), so that 

(3.1) (i?r/x)(a, b) = {RTf){a, b - a - x + x^). 





6-x 


^-^de 


(1 + 


e-x 


2ya+m-l)/2 
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Changing variables, we get 

(3.2) / {RTf){a,b) 

r(a/2) Jim 



\a ■ x' + b ■ 






(1 + 


a 


|2)(q+ 


m-l)/2 



vr{™-i)/2 



r((a + m- l)/2) 



/(l/)k-l/r Rea>0. 



Using the normahzing coefficient of the one-dimensional Riesz potential 
(cf. (O), we denote 

vi + l«l / 

X e M"", Rea>0, a 7^ 1, 3, 5, . . . ; da 



(1 + |a|2)™/2- 
Under the same assumption for a, we define 

(a, 6) G M"^. The following duality relation is an immediate conse- 
quence of Fubini's theorem: 

{R^f){a,b)ip{a,b)dadb= [ f{x){RT''^){x)dx. 



Note that 

\a ■ x' + b — Xr 



dist(x, h) 



is the Euclidean distance between the point x G and the hyperplane 
h = {y ^ M™' : a ■ y' + b = ym}- Thus, the kernel of fractional integrals 
f l3.3l) and (13.41) has a natural geometric meaning. By Fubini's theorem, 

(3.5) R^f = ijJ'^Rrf, Rt"-^ =Rt ^^AV, 
where 

V'<,(a) = (l + |a|2)(l-")/^ 

/" denotes the one-dimensional Riesz potential of order a in the b- 
variable, and 

(3.6) {Rt ^){x) = / ip{a, Xjn — o, ■ x') da. 
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We call {Rt 'p){x) the dual transversal Radon transform of (p. By 
Fubini's theorem, 

(3.7) / {RTf){a,b)<p{a,b)dadb= f f{x) {Rt <p){x) dx 

provided that at least one of these integrals is finite when / and ip are 
replaced by |/| and \ip\, respectively. Indeed, 



l.h.s = I da dx' / (p{a,b) f{x' ,a ■ x' + b) db 

f(x)dx / Lp{a, Xm — ci ■ x') da = r.h.s. 



For sufficiently good / and if, owing to the limit property of Riesz 
potentials, we can pass to the limit in (13. 5p as a ^ 0. This gives 

lim R^f = Rrf, lim Rr'^cp =Rt 
{RTf){a, b) = v/l + |a|2 (i?T/)(a, &), <^(a, b) = ^/l + H^yp(a, b). 



Theorem 3.1. Let ^ = Rrf, f G Lp{ 

m 

I <p < 

a + rr. 

Then for almost all x G 



m 

l<p< , iteQ;>0, 0:7^1^3^5^ 

a + m — 1 



(3.8) {RT'"^)ix) = (27r)™-i (r+"-V)(x). 

Moreover, for 1 < p < m/{m — 1) and almost all x G 



(3.9) {Rt ^){x) = (27r)™-i il"'-'f)ix), ^(a, b) = ./l + WAa, b). 

Proof. Equality (13. 8p coincides with (13. 2p above. Equality (13. 9p can 
be formally derived from (13. Sp by passing to the limit as a — 0. A 
rigorous proof for arbitrary / G L^(R™), 1 < p < m/{m — 1), is the 
following. Note that for such p the Riesz potential f){x) is finite 
for almost all x G M™. 

Owing to (13. ip . it suffices to prove (13.90 at s = 0. Thus, we have to 
show that 

{RTf){a,0)da _ TT^^-'y' f m 

(i+|ap)(m-i)/2 r((m-i)/2) A™ \y\ ^' 

Changing the order of integration and passing to polar coordinates, we 
obtain 

I dy' I , .„,^L_1w, dr. 



A™., (l+^2)(m-l)/2 
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Then, changing variables, we get 

m\ymr-'Ay)dy, j{y) = \ (yi + [a%nm-^)i^ - 

If m = 2, then 

1 



+ 



L(2/i+ (1-1/1)2)1/2 (^2 + (_l.^^)2)l/2j • 

If m > 2, then J(|/) can be evaluated using (12.11) and |GRy[ formula 
3.238(3)] as follows. Set A = ^ = y' /W\- Then 



J{y) 



da 



2|l/'|™-^ is—^ (A2 + {a ■ cu)2)(™-i)/2 



1 (l_t2)(™-4)/2^^ 1)/2 





2 


—m 




y 





\yT'^ Jo (A2 + t2)(m-i)/2 r((m-i)/2) 

This gives what we need. □ 

4. Inversion of the Transversal Radon Transform 

By Theorem 13.11 inversion of the transversal Radon transform Rt 
reduces to inversion of Riesz potentials of L^-f unctions on M™. The last 
topic was studied in numerous publications; see [SKMi IRulj ISam3j and 
references therein. The corresponding results are adapted in [RuT] for 
the L^-theory of the fc-plane Radon transforms on M". 

Theorem 13.11 and the second equahty in (13.51) suggest the following 
approaches to reconstruction of / from if = Rxf ■ 

The first approach: 

Invert the m-dimensional Riesz potential y in (13.81) to get 

(4.1) / = (27r) ^-"^ ©"+"^-1 i?V<^, 

where D"+''"-i = (_A)("+™'"i)/2 jg ^j^g Riesz fractional derivative. Here 
a can be chosen as we please. 

The second approach: 

Set (formally) a = 1 — m in (13. 8p to get 

(4.2) / = (271)1-"^ i^V-'^y?. 
* 

Of course, RT^~'^ip and all Riesz derivatives should be properly in- 
terpreted. 

The described method is known in the literature as the method of 
Riesz potentials; cf. |Ru7l Section 4], |Rou] . Below we provide formulas 
(14.11) and (14. 2 p with precise meaning. 
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4.1. The first approach. 

4.1.1. Hypersingular integrals. We assume a = 0, when (13. 9p yields 



/ = (27r)i-'"D™'i Rt ^, ^(a,6) = VTTH^(^r/)(a,6). 
Consider finite differences 

(4.3) {Alg){x) = J2(^)i-^y9{x-jy), 

(4.4) (a5^)(x) = J2 H i-^yaix - Vjy), 

and the relevant normalizing constants 

(4.5) dm/{m — l)= / \2m-i (^1 ^^^^ coordinate of y), 

\y\ 

oo 

2'^~"^7r^/'^ /" fl — e~*)^ 

(4.6) dmfc(m-l) = — / ^ , ,,,,, dt. 



The parameters k and £ in these formulas will be chosen as follows. 
We assume k to be any integer greater than (m — 1)/2. Furthermore, 
we set £ = m — 1 if m is even, and fix any £ > m — 1 if m is odd. 
Integrals (14.51) and (14.61) can be explicitly evaluated and the following 
statement holds; see [Rul, pp. 238, 239] and [SKM, Section 26] for J"/ 
with arbitrary a. 

Theorem 4.1. Let g = r'-^f, f G Lp{W^), 1 < p < m/{m - 1). 
Then 

•'^^ d Jm -1) l7/|2m-l 

(im,l\lll' J- J jRm \y\ 

1 [ i^Mx) 



r [A^gjix^ 

/ U,|2m-1 

JM'" \y\ 



dm,k{m - 1) 

where L™ = lim T . This limit exists in the L^-norm and in the a.e. 
sense. For / G Co fl L^, it exists in the snp-norm. 

Theorem 14.11 and (13. 9p give the following result. 

Theorem 4.2. Let f G Lp(M™), 1 < p < m/{m - 1). We set 

(4.7) g = i2ny~^ Rrip, ^(a, 6) = v/T+H^ (i?r/)(a, 6). 
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Then f can be reconstructed by the formulas 

dm,e{m-l) JiRm \y\ 



or 



1 /- (A^^?)(x) 



(4.9) /(x)=- / , r' ' dy, 

where convergence of hypersingular integrals is understood as in Theo- 
rem 



4.1 



Remark 4.3. Theorem 14.21 includes continuous functions, the behavior 
of which at infinity is specified accordingly. We denote 

CsiUn = {fe C{Rn ■■ fix) = Oi\x\-')}, 6 > 0. 

US > m — 1, then there exists p such that m/6 < p < m/(m — 1). Hence 
inversion formulas (14.81) and (14.91) are applicable to any / G Cs{M."^) 
with 6 > m — 1. 

4.1.2. Powers of "minus- Laplacian" . Another series of inversion for- 
mulas can be obtained using integer powers of the operator —A, where 
A is the Laplace operator on M™. 

Definition 4.4. For A G (0, 1), let Lip^°'=(M'") be the space of functions 
/ on having the following property: for each bounded domain fl in 
M™, there is a constant A > such that 

\f{x) - f{y)\ < A\x -y\^ '^x,yeQ (the closure of fi). 

We denote 

C|(R") = {/ : feCsiM."") n Lipf (R") for some A G (0, 1)}. 
Theorem 4.5. Let 

(7 = (27r)i-"^ Rt^, ^{a,b) = ^/uW{RTf){a,b). 

(i) Ifm IS odd, 5>m-l, and f e Q(M'"), then 

(4.10) /(x) = (-A)("-i)/2^(x). 

(ii) If m is even, 5 > m — 1, and f G Cs{W^), then 

(4.11) /M^A/ (-AfV^-M.) - (-Ar/---^ 

c"m jr™ \y\ 

where = lim Jj^|^^ uniformly in x & M'*". 

Moreover, if m > 4, / G C|(]R"'), 6 > m — 1, then 

(4.12) f{x) = -{-AY"^-'y\l'Ag){x). 
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All derivatives in ( [y^.iO| j-( f^.i^ exist in the classical sense. 

Proof. These statements are consequences of known facts for potentials 
and singular integrals. In our case g = I"^~^f; see fl3.9p . 

(i) To "localize" the problem, let x G Br = {x : \x\ < R} and 
choose x{x) e C~(M") so that 

< x{x) < 1, x{x) = if |x| < i? + 1, and x{x) = 1 if |a;| > + 2. 

We have f = h + h = xf, ^ = (1 - x)f, 

..Jo if|x|<i?+l, , fix) ii\x\< R+1, 
J^y^)-]^f(x) if|x|>i?+2, ^2la;j-|Q if|x|>i?+2. 

Let g = gi + g2, gi = /"^"Vi, 92 = I'^-'f2. Then g^ G C^{Br), and 
for all multi-indices 7, 

d'^gi{x) = — — — / fi{y)d^\x -y\'^dy. 
'jm{m - 1) J 

\y\>R+l 

In particular, for m odd, we get (— A)^"^~^^/^(7i(x) = 0. The function 
g2 belongs at least to C"^~'^{Br), and differentiation is possible under 
the sign of integration; see, e.g., [Vll Section 1(6)]. 

Hence, for m odd, {—AY"^~^^/'^g2 = 1^/2 (the Newtonian potential 
over a bounded domain), and (i) follows by Theorem 11.6.3 from [Mit 
p. 231]. 

(ii) Consider the case m even. By reasoning from above, 

and 04.111) holds owing to Remark 14.31 If m > 4 then, as in (i), we 
have -Ag = I""-^ f . Hence -I^Ag = I"'-'^f, and (KWf follows. □ 

4.1.3. Wavelet transforms. Riesz fractional derivatives 3°' can be rep- 
resented using continuous wavelet transforms; see |RuH Section 17]. 
We apply this fact in the context of the Radon inversion procedure. 

Definition 4.6. A measure /z on M*", m > 2, is said to be radial if 

g{jx)dfi{x) = / g{x)dfi{x) 

for all rotations 7 G SO{n) and all /x-integrable functions g. 
Let /i be a radial finite Borel measure on M"^ and let 

{Wg){x,t) = {g * fit){x) = / g{x -ty)dfi{y) 
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be the convolution of g with the dilated version of /x. If /i has a cer- 
tain number of vanishing moments and obeys some regularity condi- 
tions, then {Wg){x,t) is called the continuous wavelet transform of g, 
generated by the wavelet measure fi. In particular, if /i is absolutely 
continuous, that is, d^{y) = w{\y\) dy, then 

iWg)ix,t) = ^[ giy)wO-^^) dy. 



Formally, 



^g){x) = -ji-, I 2^^^^ dt, d,{a) = const. 



d^ia) 



that can be easily checked using the Fourier transform. Below we apply 
Theorem 17.10 from [Rulj to our case, when a = m — 1 and g has the 
form (1131). Recall that by dSSD, g = I'^'^f, f G Lp(M™). 

To state the result, we assume that /i is a radial Borel measure, 
satisfying 

(4.13) / |a;|^ (i|yu|(x) < oo for some /3 > m — 1, 

J\x\>l 

(4.14) / xUfi{x) = for IjI = 0,2,4, ... ,2[(m- l)/2]. 
We denote 

—m/2nl-m 

(4.15) 



^'"^ a„_ir(m - 1/2) 

r((l-m)/2) / \yr-^dfi{y) if m is even. 



X 



2f-i)('^+i)/2 



yl"^ log \y\ dfi{y) if m is odd. 



I ((m-l)/2)! 

Theorem 4.7. Let f G Lp(M™), l<p< m/{m - 1), 

g = {2ny-^ Rt^, ^{a,b) = y^l^ {RTf){a,b). 
If dn^m 7^ 0, then f can be reconstructed by the formula 

dfj,,m Jo t dfj^^rn J e t 

where the limit exists in the L^-norm and in the a.e. sense. For f G 
Cq n L^, it exists in the snp-norm. 

Examples of wavelet measures can be found in [Rull Section 17.4]. 
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4.2. The second approach: Ridgelet-hke transforms and the 
convolution-backprojection method. 

4.2.1. Preliminary discussion. Our aim is to give precise meaning to 
the formula 

(4.16) / = (27r)i-"^ Rr'-^'ip, ^ = Rrf, 

where / G L^(R™) and Rp^~"^ip is obtained by formal substitution 
a = 1 — m in the fractional integral in (13. 3p . Of course, we cannot 
replace a by 1 — m directly. To overcome this difficulty, we proceed as 
follows. Recall that 

I (3 ■ X ~\- h XfYl I 1 • , / 7 \ 

= — = dist(a;, h) 



is the Euclidean distance between the point x G M™ and the hyperplane 
h = {y : a-y' + h = y^}. Thus, 

* 1 /" ~ 

(i?T"<^)(a;) = — - / <^(a,6) [dist{x, h)]''-^ da db. 
7i(") it- 

The kernel of this integral operator has a singularity on the hyperplane 
h. Let us write this kernel in the form 

oo 





Here w{s) is a measurable function (it will be specified later) such that 
\ w{s)\/s°' ds < oo and c~^^ ^ 0. Changing the order of integration, 
we obtain 

(4.17) (i?>^)(x) = ^ fffMM,, 
where 

(4.18) (H/^)(a:,t) = l /■ ^(a,b) w ( ^^^^^ ] dadb. 



If w is a wavelet function with a certain number of vanishing moments 
and a suitable decay at infinity, then (I4.17P can be meaningful also for 
negative a, in particular, for a = 1 — m. In this case 

(4.19) {Rt'-'-^){x) = c„,. r ^^^^4^ dt, 



where Cm,w is a suitable constant and convergence of the integral is 
interpreted in a proper way. 
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An idea of this approach in a shghtly different form was suggested 
in |Rul[ Section 10.7] as a generahzation of Marchaud's method in 
fractional calculus. 

Note that / can be reconstructed directly from {WRTf){x,t) if this 
composition represents an approximate identity. This observation is 
a core of the convolution-backprojection method in tomography [Nal 
IRu4j and agrees with the method of approximative inverse operators 
in fractional calculus; cf . |NSl ISam3j . 

Our first task in this section is to show that w can be chosen so that 

(4.20) ]mi{WRTf){x,t)=c^f{x), 

where 7^ is a constant depending on w and the limit is understood 
in the L^-norm and in the almost everywhere sense. Then we replace 
the passage to the limit in fl4.20l) by integration against the dilation- 
invariant measure dt/t and get 

(4.21) r ^^^^£l}!^dt = clf{x), 4 = const 7^0. 

* 

Formula fl4.2ip agrees with equality fl4.16p . in which Rrp^~^Lp is repre- 
sented by (I4A9D . 

The convolution (I4.18P has the same nature as the ridgelet transform 
in |Cll IC2] . A detailed exposition of this topic for totally geodesic 
Radon transforms can be found in |Ru8] . 

4.2.2. Justification of U.20\ ). 

Lemma 4.8. Let f G Lp{W^), 1 < p < m/{m — 1), m > 2. Suppose 
that 

(A99^ ( ^ - S 0{s''-') if s < 1, 

w[s) I o(s-™/p'-^2) ifs>l, l/p+l/p' = l, 
for sufficiently small £i,£2 > and let 

(4.23) kiy) = ^ t'\\y\' - s^Y-'^Ms) ds, 

\y\ Jo 

kt{y) = t'^'kiy/t), t > 0. Then 

(4.24) {WRTf){x,t) = {f*kt){x) 
for almost all x. 

Proof. Changing variables and using Fubini's theorem, we obtain 
{WRTf){x,t) = [ f{x-ty)k{y)dy, 
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k{y) = I w\ '^"l^l^" 1 da. 



m — 1 



1 + |a|2 



By (123!), 



%)= / _^wi\y ■ e\) d9 = ^ I wi\ye\)d9, 



and ([2IID yields 

k{y) = am-2 I {l-ty'--''>^M\y\t)dt. 



This gives f l4.23p - fl4.24p . To complete the proof we must justify appli- 
cation of Fubini's theorem. It suffices to show that {WRTf){x,t) is 
finite a.e. when / and w are replaced by their absolute values. We 
split the integral 



da db 



(4.25) I{x,t) = ^ f {RT\ma,b) ^f^^^\ 

into two pieces Ii + I2, where in Ii we integrate over the set {(a, fe) : 
dist{x, h)/t < 1} and in I2 over the set {{a,b) : dist{x, h)/t > 1}. 
Owing to fl4:22|) and fl33D . 

/,<cti-^-'-(i?>i?T|/|)(x), 2 = 1,2; 
Pi = El, (32 = 1 — m/p' — £2, c = const. 



By Theorem 13. U are dominated by multiples of Riesz potentials, 
namely, 

(4.26) h < ct^-^'-"^(J'^'+"^-Vl)(2;)- 

The latter are finite for almost all x (see, e.g., ^Stlt Chapter V, Section 
1.2]) because, for sufficiently small Si and £2 > we have Q < (3i+m — 
1 < m/p. This completes the proof. □ 

Owing to Lemma 14. 8[ the standard machinary of approximation to 
the identity (see, |Stlt Chapter III, Section 2.2])) implies the following 
statement. 

Theorem 4.9. Let 1 < p < m/{m — 1), m > 2. Suppose that w{s) 
satisfies l!i4-22 ) and the corresponding kernel 

/-l^l ,2 



%) = ^ / {\y\' - s^Y-'^Ms) ds 
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has a radial decreasing majorant in L^(R™). If f & L'^iW^), then 
(4.27) lim-^ / {RTf){a,b)w r'f^llZlA dadb = jf{x) 



(4.28) 7= / k{y)dy, 

where the limit exists in the U'-norm and in the a.e. sense. If, more- 
over, f G Co(ffi™) nL^iW"), then the convergence in ^-STl ) is uniform 



on R"". 



To give an example of a function w satisfying Theorem l4.9[ we invoke 
Riemann-Liouville fractional integrals (12.51) and write k{y) in the form 

k{y) = 7r("*-i)/V2-™(4';^-')/'[s-i/2^(v^)])(r2), r = \y\. 

Then 

r poo 

7= / k{y)dy = c / \{t)dt, 

JR™ Jo 

m-1/2 

A(t) = (4r)/^[s-V^^(v/i)])(t), c=^^. 
Example 4.10. (cf. [SS], p. 338). Let 



r(i + a) it + iY+^' 

Both functions are integrable on (0, oo) for £ > a and 

/•OO 

(4.29) / Xa/it)dt = T{i- a). 

Jo 

We set a = (m — l)/2. Then for any i > (m — l)/2 the function 



Wis] = SK 



{m-l)/2,e\ 

can be used in Theorem 14. 9[ The corresponding constant 7 in (14.281) 
will be computed as 

r,,,, r, 7r'"-i/2r(£- (m- l)/2) 

7 = c / X{t)dt = c / Xi^m-i)/2A't)dt- ^ ^ ' 



"^^^-^^ - r(m/2) 
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4.2.3. Justification of Let us represent the corresponding trun- 



cated integral as an approximate identity. 

Lemma 4.11. Let f E L^iW^), 1 <p < m/(m — 1), m > 2. Suppose 
that w satisfies and let 

aA.y) = e^"'s(y/e), £ > O. Then 

(4.31) f(lMli^,,.U*.M^) 

for almost all x. 

Proof. The result follows from (14.241) if we change variables and apply 
Fubini's theorem. To justify application of Fubini's theorem, we note 
that the left-hand side of (I4.3ip is finite a.e. when / and w are replaced 
by 1/1 and respectively. Indeed, let 

U^) = j'^^-^dt^ 

where I{x,t) is defined by K2^ . By (|4:26|) . 

' dt 



fi3i+m ' 



Ux)<cY^{I^^^--'\f\){x) / 

i=l ^'^ 

which is finite a.e. because (3i + m> 1. 

It remains to perform simple calculations. We have 

Set k(y) = kMy\). Then, by 623, 



t \ i J 

where g has the form (14.301) . □ 

Lemma [4.111 implies the following inversion result. 

Theorem 4.12. Let f G Lp(M™), 1 < p < m/{m - 1), m > 2. 
Suppose that w{s) satisfies ^4-2^ and the corresponding kernel 

9{y) = ^^^^ /"Wp - s'Y-^yMs) ds 



(m-l)|y|-7o 
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has a radial decreasing majorant in L^(M™). // 

{Wif){x,t) = ^ [ <^{a,b)w ( 1 dadb, 

^ jR™ \ tv 1 + lo-r 

where (p = Rrf , then 

(4,32) ..^^jr-OMM,,., 

the limit being understood in the U'-norm and in the a.e. sense. If, 
moreover, f G Co(]R™') fl L^(R™), then the convergence in ^4-32^ is 
uniform on . 

To give an example of a function w satisfying Theorem 14 .121 we write 
g{y) in the form 

{m-l)/2 

9iy) = {lt^'^''[s-''M^mr\ r = \y\. 

Then 



(4.33) 71= /" g{y)dy = c [ Xi{t)dt, c 

^R™ Jo 



2r(m/2)' 



(4.34) X,{t)=t~\l!-^'^/'[s-'/MV^)]m- 
Example 4.13. Let 

w(s) = SK(„+i)/2/(s^), ^>(m-l)/2; 
see Example 14.101 Then 



Ai(t) = t ^A(m,+l)/2/(i) 



i(m + 1) 

Hence, by ( Km . 



r(i + (m + i)/2) {t + iy+^ 

\rn-l)/2A^) (GL1(0,Oo)). 



= c I Ai(t)rft = ^^^^^^ A(^_i)/2,,(t)cit 







7r™-i/2r(£- (m-l)/2) 
i(m + l)r(m/2) ' 
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Remark 4.14. The assumption g G L^(R™) in Theorem 14 . 1 21 can be ful- 
filled only if the generating function w has vanishing moments. Hence 
{W(p){x, t) is a wavelet-like transform (or a ridgelet-like transform) and 
fl4.32p has the same nature as the classical Calderon identity; cf. |F JWj . 
Furthermore, by Lemma [2.21 for any function w, satisfying 



=2i 



w{s)ds = Vj = 0,2,...,2[(m- l)/2]. 



7i = c < 



(is < oo for some P > m — 1, 

the corresponding function (14. 34^ has a decreasing integrable majorant. 
In this case the constant 71 in fl4.33p has the form 

2T ( — - — J J s^^^w{s) ds if m is even, 

4/ ]_\(m+l)/2 roc 

/ s"^~^w(s) logs ds if m is odd. 

U(m-l)/2)!yo 

Such functions w{s), having fast decay as s and 00, can be con- 
structed using Example 1.6 from |Ru3j . 

A similar remark addresses to the assumption k G L^{W^) in Theo- 
rem U]9l 

5. Transversal Radon transforms on the 
Semyanistyi-Lizorkin spaces 

Our aim in this section is to bring light to Strichartz's inversion 
formulas (11.31) in the general context of the transversal Radon transform 
and to determine spaces, say, X and Y, of smooth functions such that 
Rt acts from X onto Y as an isomorphism. 

5.1. Isomorphism of Rt- To start with, we observe that if / G 
//-"^(R™) is nonnegative, then ip = Rxf is integrable on R'" only if 
/ = 0. This is obvious from the equality 



(/9(a, h) dadb- 



: / da db f{x',a-x' + b)dx'=\\f\\i / da. 



Thus, the relation RtJ G L^(R™) is possible only if / is sign-changing. 
It means that embedding i?r[S'(M™')] C S'(]R'") does not hold, although, 
i?T[^(M™)] C C~(M'"). 

A close situation is known in the theory of operators of the potential 
type, where Semyanistyi-Lizorkin spaces of "good" functions come into 
play and lead to deep results; see |RuH ISani3l ISKM] . These spaces were 
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introduced by Semyanistyi ^Selj and essentially generalized by Lizorkin 
[Lizlj - [Liz3] and Samko |Samlj - [Sam3j . They incorporate Schwartz 
functions, vanishing on a given set, and their Fourier images with zero 
moments. 

Given a function f{x) = f{x',Xm) on R™, let 



f{x)e'''-ydx 



be the Fourier transform of /. We denote by jFi and JF2 the similar 
transforms in the x'-variable and the x^-variable, respectively. Let 

= \E'(]R™) be the subspace of functions ip G S'(M™') vanishing with 
all derivatives d^ip, /c G Z+, on the hyperplane Xm = 0. We denote 
by $ = $(M™)^= the Fourier image of ^(R'"). The spaces 

\1'(R™) and $(R™') have the same meaning. 

The following auxiliary statements are immediate consequences of 
preceding definitions. For convenience of the reader we present them 
with proofs. 



Proposition 5.1. 

(i) A Schwartz function belongs to $( 



if and only if 



(5.1) 







yk G Z: 



Wx' G 



a'/n— 1 



(ii) J-'i is an automorphism o/\l'(R™). 

(iii) J-'2 acts as an isomorphism from \1'(R™') onto 

(iv) The map : \1'(R™) —>■ "^{W^) is continuous for every j G 



Proof, (i) Let G $, that is, 
k G Z+, 



J-'i/) for some ip & Then for any 



1 X , Xr 



'ZXr 



Owing to injectivity of J-'i^, the result follows. 

(ii) Let ip & tpi = J-'iip. Then for any k G Z+, 

{dtijM,0) = [ {dt^){x',0)e'''-y'dx' = 0, i.e., ^1 G ^. 
Conversely, we have = J-'iip2, where 

My',ym) = iJ'r'i'W,ym) = (27r)l''"(J-i^)(-t/',t/„) G VI/. 

(iii) If ^ G ^, then J'aV' = ^^0, ^0 = T'^T^'i) = J^fV e ^ by (ii). 
Hence J-'2ip G $. Conversely, let G $, that is, = J^ip, ^/^ G vf. Then 
= J^2^ii^ with J^ii: G ^. 
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(iv) Let e ^, that is, tpiy) = {J^-^ip){y) = (27r)-"(J^0)(-?/), e 
Then = (27r)-™(-l)ljl {d^J^(j)){-y), and therefore, by (i). 



X / {iXmY"^ (j){x , c^^^m = 0. 



The continuity of the map : \1/ — \1/ follows from its continuity in 
the topology of □ 

Theorem 5.2. The transversal Radon transform Rt acts as an iso- 
morphism from $(R™) onto $(]R''"). 

The proof of Theorem 15.21 relies on the following two lemmas. 

Lemma 5.3. For any f G L^{W^), 

(5.2) J-2[(/?T/)(a,-)](0 = V(a,0 G . 

This formula can be easily obtained by direct calculation, using Fu- 
bini's theorem. Introducing a "mixing" map 

(5.3) (An)(a,0=M(-ae,0, 
we formally have 

(5.4) Rrf = T:,^kTf. 

Lemma 5.4. The map A acts as an isomorphism from \1/(M™) onto 

Proof STEP 1. Let u{y) = u{y',ym) G ^(M™). For every p,q e Z+ 
there is a constant Cp,, such that 

(5.5) \uiy)\ < Cp,, \yj'P (1 + lyH"" Wy e M™; 

cf. |Sam3l Ch. 2, Section 4]. If \ym\ > 1 this inequality is obvious and 
holds for every Schwartz function. In the case \ym\ < 1 the estimate can 
be obtained by making use of Taylor's formula with integral remainder. 
We first show that 

(5.6) v{a,0 = u{-a^,Oe^i^n- 

Clearly, f gC°°(R™') and a simple calculation shows that every deriv- 
ative d^v, j G Z™, can be represented as a finite sum of expressions of 
the form 

Qi being certain polynomials. Since differentiation preserves the space 
^'(M™') (see Proposition 15.11 (iv)), f/i(a,0) = 0, and multiplication by 
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a polynomial preserves \E'(M™), it suffices to show that the function 
= is rapidly decreasing. In other words, we have 

to show that for any k G Z+ there is a constant Ck such that 

(5.7) (i + |a|2 + eYk(«,OI <c. V(a,Oet™. 

This would imply t'GS'(]R™). We have 

\2\i tc2\k-i 



l.h.s. = J2('')vk,i, V, 



i, = {i + \a\'y{er-'\v-M,o\- 



Owing to (15.51) with p = k and q = 2k — i, we obtain 

Vk,i S Ck,2k-i . , ..o , I ^10 ) — . . \ ^19 ) < Ck,2k- 



k 

This gives (15.71) with = ^ {^)ck,2k-i, which implies (15.61) . 

STEP 2. Now our task is to prove that every function v G ^'(M™) 
has the form f (a, ^) =m(— a^, ^) for some function u G \I/(M"'). We set 

(5.8) u{y) = u{y',y^) = v{-y'/ym.,yni), Vm 7^ 0; u{y',0) = 0. 

Let us show that this function belongs to \E'(M™). As above, for every 
p,q E there is a constant Cp^g such that 

(5.9) \via, 0)1 < Cp,, leP^ (1 + lap + e^)"" V(a, G M'". 

Furthermore, for ym 7^ 0, every derivative {d^u){y',ym), j G Z", is a 
finite sum of expressions of the form 

Ui{y',ym) = Qi{y', 1/2/m) {d'v){~y'/y^,ym), i G Z"^, 

where Qi are certain polynomials. Again, since differentiation preserves 
the space \E^(]R™) and multiplication by Qi{y', 1/ym) preserves \I/(M"'), 
it suffices to show that the function Ui{y) = (d^v) {—y'/ ym, ym) is rapidly 
decreasing when \y\ ^ 00 and ?/m — > 0. 

Let us check that for any r, s G Z+ there is a constant Cr,s such that 

(5.10) y-''{l + \y\'y\u;{y)\ < c.,, Vy G M"^. 

r 

We have l.h.s. = ^ (Qf^.fc, where 

fe=0 

Ur,k = y-'' (1 + yif ilyVr' Wv)i-y'/ym,ym)\. 
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Owing to (15.91) with p = s and q = 2r — k, 



Ur k ^ 2r—k 





\y'\ 




1+y 


2 

m 




\y'\ 





r—k / 7/2 r—k 

am 



1 + W^IV' 



This gives (I5.10p with Ck = Yl (fc)cs,2r-fc, which imphes u E ^(M"^). 

k=0 

To complete the proof we note that the equahty f (a, ^) = u{—a^,C, ) = 
{Au){a,^) for the function u defined by (15. 8p is obvious. Thus the map 
A : \I/(]R™') \1/(]R''") is surjective and the inverse map 

is well-defined by 

(5.11) {A~'v){y)=v{-y'/y^,ym), ym^O; {A~'v){y',0) = 0. 

Both maps A and A~^ are continuous in the topology of the Schwartz 
space. This can be proved directly using (15.51) and (15. 9p . □ 

Proof of Theorem \5.S\ Since the maps 

jr . $(K'") ^ ^(R'"), A : ^(R") ^ ^(R"^), J^:^^ : ^(R*") ^ $(R'" 

are isomorphisms, then formula (15. 4p is well-justified on functions / G 
$(R™'), and the result follows. 

5.2. Inversion formulas. A pointwise inversion formula 

(5.12) f = J^-'A-'j^2RTf, /G<I>(R'"), 

in terms of the Fourier transforms follows from (15.41) . Below we obtain 
alternative inversion formulas, that do not contain the Fourier trans- 
form. To this end we invoke the partial Riesz potential 

(rcfv \ 1 f f{x',ym)dym 2-7rV2r(a/2) 

(5.13) {I^f){x) = ^- —, 7i(«) - 



Ilia) J ' r((l-«)/2)' 

i?e a > 0, a 7^ 1, 3, 5, ... . 
This operator is an automorphism of $(R'") |SKM] and 

(5.14) T{l^f\{y) = Ivmrmiy), / e 

The last relation extends /^/, / G $, to all a G C as an entire function 
of a. We also introduce a backprojection operator which sends 



^We do not call it the dual operator because the latter has another meaning in 
our paper. 
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functions on to functions on by the formula 

(5.15) [Rj.g)[x)=l g{a,-x' ■ a + Xm)da= {RTg){-x',Xr,i). 

Theorem 5.5. If f E $(R™), then for any complex a and (3, 

(5.16) {i:^RTl^^RTf){x) = (2vr)™-i (7^'^+'"- V)(x), x e M'". 
Proof Let g = I^Rrf- Then g G $(M™) and (O yields 

(j-2/2"[(i?T^?)(x',-)])(^) = \v\-''{WTg){x',-)m 

(5.17) = |r7r"(J-^7)(x'77,r7). 

Furthermore, denoting by J^i the Fourier transform in the first m — 1 
variables, we have 

(.F^7)(x'r/,r/) = {Tr{:F,[g{a,-)]{iimx'r^) 

= {J-dy'2[{I^RTf){a,-)]{vmx'v) 
= {J^^M-^{J^f)i-arj,v)})ix'r^) 

[ {J'f){-av,v)e'^-''^'da 



= \V\ 

= Ir^l^-^-™ / {J^f){w,r])e-''"-^' dw 

= {2nr-'\v\'-'-^iJ^2[fix',-)m. 

Hence, by ( ISTTj) . 

(^2/2"[(i?T^)(x',-)])(^) = i2nr-' \v\'-''-^-"'iHfi^'rmv), 
which implies f l5.16p . □ 

Equality (I5.16P gives a variety of pointwise inversion formulas. Let 
©2 = -^2"" be the corresponding partial Riesz fractional derivative, 
which is well defined on functions / G $(]R'") in the Fourier terms: 

(.FD-/)(y) = iy^r(.F/)(i/). 

Corollary 5.6. Let (f = Rrf, f G $(]R™). The following pointwise 
inversion formulas are contained in Ii5.16\) : 

(5.18) / = (27r)^~'"©^-^/?T</', (/? = 0, a = 1 - m); 

(5.19) / = (27r)i-'" ^tD^-V, {a = 0, (3 = 1 - m); 

(5.20) / = (27r)^-'"Df''^/'i?Tm)2""'^^V, (a = /3= (l-m)/2); 

(5.21) / = (27r)i-™D™i?T/2V, (« = -m, /? = 1). 
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The Riesz derivative D2, corresponding to \ym\ in the Fourier terms, 
can be expressed through the usual derivative dm by the formula D2 = 
iH2dm, where H2 stands for the Hilbert transform in the last variable: 

{H2ip){a,b)=p.Y.— [ ^4~^dbi, (J^/^z^) (y) = sgny„ (J^y?) (y) ; 

see, e.g., |Graf] . In particular, for m = 2n + l, n G N, we have the 
following 

Corollary 5.7. Let ^ = Rrf, f e ^(M^^+i). Then 

(5.22) / = {27rr'-{-l)-dl:^,RT^ 

(5.23) = {2nr'-i-irRTdl:^,^ 

(5.24) = (27r)-2«(-l)"92Vii?T5L+i^. 

In fact, fl5.23p and fl5.23p follow from 05.221) because Rt commutes 
with differentiation in the last variable. 

Remark 5.8. Inversion formula fl5.12p and those in Corollaries 15.61 and 
15.71 are pointwise analogues for Rt of the Strichartz's formulas fll.3p . 
In the next section we adapt these formulas for the Radon transform 
on the Heisenberg group. 

6. Connection between the Radon transforms R, Rt, Rh 
We denote 

a;= (xi, . . . , x„) = (x', Xm)eW^, 9 = {9^, . . . ,9^) = {9', 9^) G 

and let h he a hyperplane which is not parallel to the The 
following two parameterizations are available: 

(6.1) h = {x:xm = a-x' + b}, (a, 6) G M"^ = M"^"^ x M, 

(6.2) h = {x:9-x = t}, (0, t) G P"" = X M, 0^ ^ 0. 
Different sets of parameters are related by 

(6-3) e= t ^ 



(6.4) a = -f, b = ^. 

The corresponding transference operators acting on functions ip{a, b) 
and ij^{9,t) have the form 

(6.5) {Tv){0,t) = ^(-^,^ 
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(6.6) 



(r-V)(a,6) = V' 



a - 



\a\^ a/1 + |a,| 
Lemma 6.1. The following relations hold: 

(6.7) {Rf){0,t) = \e^r\rRTf){e,t), ^ o, 



(6. 



{RTf){a,b) = {l + \a\Y'^\T-'Rf){a,b). 



Ifilj{6,t) = ip{—0, —t) for all {6,t) G P™, then for the dual transforms 
M.12\) and ( TT^) we have 



(6.9) 

Furthermore, 
(6.10) 



{R*il)){x) = 2{Rt r-V)(a;). 



{Rt ^){x) = -{WT^){x). 



Proof. Equalities (ETD and (ESD follow from ffLT3ll . Equality (ES]) can 
be obtained by making use of (12. 3p . Namely, 



(i?*^)(a;) 



i){e,d-x)de 



a + Cr, 



{a + Cm) ■ X 



+ ^ 



a - 



a — 



a - 



X 



da 



a + |a|2 Vl + |a| 
2 / (T-^^lj){a,Xm-a-x')da = 2{RTr-^^lj){x) 



Equality ( 16.10p is a consequence of ( 16. 9p . It can also be obtained di- 
rectly by using (12. 4p . □ 

Connection between the Radon transforms Rt and Rh is given in 
(11.61) . Let us present this equality and a similar one for the dual trans- 
forms in the operator form. We define 

(i„,)(C.r, ^ 2/^^,((C..,o(.,0))|j^^^l^ 

, 1 r -XX T T 2^2; 

^(^(2;, T--Im iC-z)) dz, dz = ^^^|^|2)n+l/2 ' 



(6.11) 
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and set 

X = (x(i),a:(2),X2„+i), = (xi, . . . ,a:„), X(2) = (x„+i, . . . , X2„), 

a = (a(i), a(2)), a(i) = (ai, . . . , a„), a(2) = (ctn+i, . . . , a2n)- 
Lemma 6.2. Given functions f and (f on EI„, let 

(6.12) ( Q/) (x) = /(x(i)+ix(2) , X2„+i ) , ( Q(f) (a, b) = V9(2a(2)-2m(i) , 6) . 
Then 

(6.13) = RtQ. 

(6.14) QRh=RtQ. 

The duality relation 

(6.15) / /(C,r)(iH</^)(C,r)^^Ct^^= / {RHf){z,t)^{z,t)dzdt 

holds provided that at least one of these integrals is finite when f and 
(f are replaced by \ f\ and \ip\. 

Proof. Equality ( ]6.13p follows from ( ]1.4p -( fL6l) . Furthermore, by ( 13. 6p . 
{Rt Q^)ix) = / {Q(p){a,X2n+i - a ■ x')da 
V9(2a(2) - 2m(i) , X2„+i - a ■ x') 



1 + |aP)"+i/2 



, 1 , dudv 

2 ip{u + IV, X2n+1 - 7t(^ ■ ^(2) - V ■ X^i)^ 



2' (4+ |m|2 + |t;|2)"+i/2- 

Setting X2„+i = r, X(i) = ^, X(2) = rj, ^ + ir] = (, u + iv = z, we get 

V9(2;,r - ^/m(C ■ 2))(i2;. 



This gives f l6.14p . The duality relation (16. 150 follows from (13. 7p by the 
same reasoning. □ 
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7. Inversion of the Heisenberg Radon Transform 

Equalities fl6.13p and (16.141) enable us to convert inversion formulas 
for the transversal Radon transform Rt (with m = 2n + 1) into those 
for the Heisenberg Radon transform Rh- We skip routine calculations, 
which mimic the proof of Lemma [6 .21 As in (11.11) . we write (C? t) for 
the argument of / and {z,t) for the argument of Ruf- 

Let us start with Theorem 14.21 and set 

(7.1) g{C,r) = {2T:r^-{RHm,r), 
i,{z,t) = ^l + \z\^/A{RHf){z,t)- 

(7.2) {^Ug){C.r) = j2{%-^y9{C-Jz.r-jt), £ > 2n, 

(7.3) (Af,,)^?)(C,r) = X;(^V-l)^-^(C-V7^,r-V7t), k > n- 

oo 

"^"'^ = L.rW^"^^' "^"'^ ^ r(2n + l/2) j 



Theorem 7.1. Lei / G Lp(H"), 1 < p < 1 + l/2n. Then f can be 
reconstructed by the formulas 



(|2P + i2)2»+V2 



or 

1 /■ (Af.,.)9)(C. 



where = lim^^o i|2|2_,_(2>£. ^/ie limit exists in the LP-norm and in 
the a.e. sense. For f G Co(IH[") f]L^(]HI"') it exists in the snp-norm. 

The next statement follows from Theorem 17.51 

Theorem 7.2. Let f G Q(H„), 6 > 2n, and let g be defined by ( fTJ] ). 
Then 

(7.4) f^^ + ,^^r) = {-Arg{^ + zr^,T), 

where 
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Theorems 14.71 14.91 and 14.121 can be reformulated in a similar way. 
We leave them to the interested reader. 

Another series of pointwise inversion formulas for Rh can be derived 
in the framework of the corresponding Semyanistyi-Lizorkin spaces. 
Following Section 15. we define 

(7.5) $(H„) = {0 G 5(e„) : / 0(C,r)r*^rfr = O VA; G Z+, VC G C"}. 



Theorem 7.3. The Heisenberg Radon transform Rh is an automor- 
phism of the space $(EI„). 

Proof. By f l6.13p . Rh = Q~^RtQ- Hence, Theorem 17.31 is a conse- 
quence of Theorem 15.21 □ 

Theorem 7.4. A function f G $(IH[„) can be reconstructed from its 
Heisenberg Radon transform if = Rnf by the formula 

(7.6) /(e + IT,, r) = — ^ / wiy) e-^ («''''-) dy, 
where 

w{y) = U{- e^^^-+i dh, y = (i/(i), ^(2), l/2„,+i). 

Proof Replace / by Qf in fl5.12p and make use of fl6.13p . This gives 

which coincides with ( 17. 6p . □ 
Inversion formulas below do not contain the Fourier transform. 

Theorem 7.5. A function f G $(IH[„) can be reconstructed from ip = 
Rnf by the following formulas: 

(7.7) /(C,r) = (-l)"(4vr)-2"a^(i?^^)(C,r) 

(7.8) = (-l)"(47r)-2"9^(i?^,9,2»(C,r) 

(7.9) = {-iriAn)-'-d-{RHd-^){C,T). 

Proof. Since Rh commutes with differentiation in the last variable, it 
suffices to check the first formula. We replace / by Qf in fl5.22p and 
make use of fl6.13p . This gives 

/ = (-1)" (2vr)-2" Q-'di:^,RTQ^. 
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(Qv')(a, -x' ■ a + X2n+i) da 
= / v^(2a(2) - 2m(i), -x' ■ a + X2„+i) (ia 

= 4""/ + if , -^(f ■ - M ■ X(2)) + X2n+l) C^wdf . 

Setting C = ^ + z?7, 2; = M + it>, we have 



This gives the resuh. 



4-" / v^(^,r-^/m(C-^))rf^ = 4-"(i?^^)(C,r) 



□ 



Theorems 17.41 and 17.51 give precise meaning to Strichartz's formulas 
in ([L3D. 

8. Mixed-norm estimates for Rt and Rh 
For 1 < g, r < 00, we define the following spaces with mixed norm: 



i^(^,t)rcit 



-1 (j/r 



1/9 



< 00}; 



= {i^ia,b) : 
= {9{z,t): 

Wdi II g,r 



|V(a,6)r'(i6 



q/r 



1/g 



(ia I < 00}; 



\g{z,t)\'-dt 



q/r 



1/g 



dz I < 00}. 



Theorem 8.1. For m > 2 an a priori inequality 

\\Rf 1^ \\q,r — Cp^q^r II / lip 

holds if and only if I < p < m/{m — 1), q < p' {1/p + l/p' = 1) and 
1/r = 1 — m/p' . 
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This well-known statement, which is due to Oberlin and Stein |0S] 
and also Strichartz [Strlj , implies the following result for the transversal 
Radon transform. 



Theorem 8.2. For m > 2 an a priori inequality 

||i?T/;K 

holds if and only if 
(8.1) 1 <p < m/{m - 1) 



\q,r — "--p,g,r II J ||p 



q=p 



1/r = 1 — m/p . 



Proof. Note that, unlike Theorem 18. now we have a strict equality 
for q. Using (16. 8p . (16. 6p . (12. 4p . and taking into account that 



{Rf){-e, -t) = (Rfrnt), 



we obtain 



q,r 



1 + la 



2N-1/2 



m 



a — Cr 



JR™-! '-Jr ^ V t 



1 + |aP 

q/r 



db 



dt 



1 + |a|2- 
(l + |a|2)(^-^')^/2-da 



q/r 



da 



\{Rf){eM dt 



q/r 



de. 



Now we set q =p\ 1/r = 1 — m/p', and apply Theorem 18.11 This gives 



iRrf;'. 



I q,r 



2-^/''\\Rf;F"'\\g^r < 



-1/9, 



•^p,q,r II J IIP) 



as desired. 

The necessity of (18. ip can be proved using the homogeneity argu- 
ment, as in [5tr2j . □ 

The mixed-norm estimate for the Heisenberg Radon transform now 
follows by Lemma 16. 2[ 



Theorem 8.3. For n > 1, an a priori inequality 

\\RHf',^n\\q,r < Cp,q,r' || / ||p 

holds if and only if 

(8.2) l<j9<l + l/2n, q=p, 1/r = 1/p - 2n/p' . 
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Proof. We have 

(set u = 2a(^), 



\{RHf){u + tv,t)\'dt 



1 



\{RTQf){a,b)\'db 



dudv 
b) 



q/r 



da. 



iP(H„)- 



Hence, by Theorem I8.2[ 

IIP f ■ THT II <^ AJ^/l r 

This gives the result. □ 

Theorems 18.21 and 18.31 generahze Corollary 4.2 in jStr2j and extend 
it to the full range of p; cf. Theorem 11.31 and discussion after it. 
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